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ABSTRACT 



Context. Stellar flares affect all atmospheric layers from the photosphere over chromosphere and transition region up into the corona. 
Simultaneous observations in different spectral bands allow to obtain a comprehensive picture of the environmental conditions and 
the physical processes going on during different phases of the flare. 

Aims. We investigate the properties of the coronal plasma during a giant flare on the active M dwarf CN Leo observed simultaneously 
with the UVES spectrograph at the VLT and XMM-Newton. 

Methods. From the X-ray data, we analyze the temporal evolution of the coronal temperature and emission measure, and investigate 
variations in electron density and coronal abundances during the flare. Optical Fe xiii line emission traces the cooler quiescent corona. 
Results. Although of rather short duration (exponential decay time t^c < 5 minutes), the X-ray flux at flare peak exceeds the quiescent 
level by a factor of a; 100. The electron density averaged over the whole flare is greater than 5 • 10" cm"'. The flare plasma shows 
an enhancement of iron by a factor of ^2 during the rise and peak phase of the flare. We derive a size of < 9000 km for the flaring 
structure from the evolution of the the emitting plasma during flare rise, peak, and decay. 

Conclusions. The characteristics of the flare plasma suggest that the flare originates from a compact arcade instead of a single loop. 
The combined results from X-ray and optical data further confine the plasma properties and the geometry of the flaring structure in 
different atmospheric layers. 
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1. Introduction 

Flares are among the most prominent signs of stellar activity. 
They can be observed over the entire electromagnetic spectrum 
from X-ray to radio wavelengths, demonstrating that all layers of 
the stellar atmosphere are affected, with a wide range of plasma 
temperatures and densities involved. A flare is a highly dynamic 
event, revealing the complex behavior of the stellar atmosphere 
with rapidly changing physical conditions, i. e. its interplay with 
varying magnetic fields in time and its response to the sudden re- 
lease of large amount s of energy. In the commonly accepted pic - 
ture of a stellar flare (iHaisch et al.lll99lHPriest & Forbesll2002l) . 
a magnetic reconnection event in the corona drives the accelera- 
tion of particles that radiate in nonthermal hard X-rays and radio 
gyrosynchrotron emission as they spiral down along the mag- 
netic field lines. Lower atmospheric layers are strongly heated 
by the impact of these particles and immediately start cooling 
by radiation in the optical and UV continuum as well as in chro- 
mospheric and transition region emission lines. Chromospheric 
evaporation then brings "fresh material", emitting soft thermal 
X-ray emission, into the corona. The strong X-ray and UV 
radiation field finally in duces further chromospheric emission 
(lHawlev&Fisheilll992h . 



* Based on observations collected at the European Southern 
Observatory, Paranal, Chile, 077. D-OOl 1(A) and on observations ob- 
tained with XMM-Newton, an ESA science mission with instruments 
and contributions directly funded by ESA Member States and NASA. 



Flares on the Sun and on stars show a wide variety of am- 
plitudes, from the smallest micro/nanoflares to giant flares with 
luminosity increases by orders of magnitude, and timescales 
ranging from a few seconds up to several days. Depending on 
energy budget, decay time and shape of the lightcurve, impul- 
sive and gradual flares can be discerned. On the Sun, the for- 
mer have been associated with compact emission regions, i.e., 
single loops, while the latter typically involve a series of erup- 
tion s in a whole arcade o f loops (so-called two-ribbon flares, 
see iPallavicini et al] |l977!). For stellar flares, the loop geome- 
try can usually not be observed directly, but scaling laws based 
on the hydrostatic case, relating the loop temperature and pres- 
sure with the size of the loo p, as derived and tested for the 
quiet Sun (Rosneretal. 1978) can be adapted to stellar flares 
to give an esti mate of the dimensions of the involved coro- 
nal structures dAschwanden et all l2008h . Additionally, hydro- 
dynamic loop modeling approaches allow to assess the loop 
length from the analysis of the decay of the fl are in X-rays and 
to estimate the am ount of additional heating (ISerio et al.lll99ll: 
iReale et aDll997h . This approach has re cently been extended to 
an analysis of the rise phase of the flare (lRealell20"07h . Individual 
flare events can thus be characterized with respect to to the phys- 
ical properties of the flaring structure from observational quan- 
tities also available for stellar events. The greater the spectral 
coverage and the better the spectral and temporal resolution, the 
more information can be gained and the better are the constraints 
on the physical properties of the flare plasma and the geometry 
of the flaring active region. Simultaneous observations in multi- 
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wavelength bands allow to study the flare throughout all atmo- 
spheric layers from the photosphere into the corona. 

Strong flares dominate the overall flux level in all wave- 
length bands; the X-ray luminosity can increase by more than 
two ord ers of magnitude during the strongest events (e. g. 
Schmitt & Favatal 119991: iFavata & Schmittl 119991: iFavata et alj 
2000h . which allows the strongest events even on relatively faint 
targets to be investigated in great detail. For outstanding flares 
like the long-lasting giant flare observed on Proxima Cen with 
XMM-Newton, the temporal evolution of plasma temperatures, 
densities and abu ndances during different flare phases could be 
observed directly dGiidel et al.l2002ll2004l) : detailed and specific 
hydrodynamic loop modeUng revealed a complex flare geome- 
try with two loop p opulations and differing heating mechanisms 
dReale et al.ll2004l) . 

Even for much smaller flares, low- and intermediate- 
resolution X-ray spectroscopy makes it possible to follow the 
evolution of plasma temperature and emission measure during 
flare rise and decay. The temperature of the flare plasma is typi- 
cally observed to peak before the maximum emission measure is 
reached, verifying the overall framework of initial heating, evap- 
oration, and conductive and radi ative coohng as predicted by 
hydrodynamic loop models (e. g. lRealell2007h High-resolution 
spectroscopy is needed to directly measure coronal densities dur- 
ing flares and thus to obtain the weights of changes in plasma 
density or volume. In general, coronal densities deduced from 
the ratio of the forbidden and intercombination lines in He-like 
triplets, especially from O vn and Ne ix at formation tempera- 
tures of 2 MK and 3.5 MK, respectively, during smaller flares are 
consistent with the quiescent val ue within the errors but tend to 
be higher (see e. g. iMitra-Kraev & Ness 2006): unambiguously 
increased coronal densities have so far only been measure d dur- 
ing the giant flare on Proxima Centauri (iGiidel et alj|2002l) . 

Another aspect of the evaporation scenario is the observa- 
tion of abundance changes during flares, leading directly to the 
question to w hat extent flares influence the ov erall coronal abun- 
dance pattern. iNordon & Behad (I2007ll2008h systematically an- 
alyzed abundance trends in a sample of flares observed with 
Chandra and XMM-Newton and found that the most active stars 
with a pronounced inverse FIP effect during quiescence com- 
pared to solar photospheric values showing a shift to a FIP effect 
compared to the corresponding quiescent values, and vice versa 
for the less active stars showing the solar-like FIP effect during 
quiescence. This confirms the existence of the inverse FIP ef- 
fect and of strong discontinuities in the abundance patterns of 
corona and chromosphere/photosphere also for stars where pho- 
tospheric abundances are difficult to measure, however, the ques- 
tion what causes this bias is left unanswered . 

We present the observation of a giant flare on the active 
M dwarf CN Leo with XMM-Newton and the UVES spec- 
trograph at the VLT. An extremely short impulsive outburst 
at the beginn i ng of the flare has already been discussed by 
ISchmitt et all (l2008h . This paper is the concluding third paper 
of a ser ies analyzing the whole flare event in greater detail. 
Paper I (iFuhrmeister et al.l l2008h reports the observational re- 
sults on the chromospheric emission as diagnosed by the UVES 
spectra, while Paper II (Fuhrmeister et al. 2010) describe model 
chromospheres and synthetic spectra obtained with the atmo- 
spheric code PHOENIX. Here we focus on the temporal evo- 
lution of the coronal plasma temperature, emission measure and 
abundances determined by the X-ray data and their implications 
on the geometry of the flaring structure, on the physical condi- 
tions dominating the flare environment, and on the abundance 
patterns predominating in different atmospheric layers. 



2. The active M dwarf CN Leo 

CN Leo (GJ 406, Wolf 359) is, after flie a Centauri system and 
Barnard's star, the fifth-nearest star to the Sun at a distance of 
2.39 pc. Despite its proximity, it is a faint object with a V band 
magnitude of 13.54; recent spectral classificati on ranks it as an 
intermediate M dwarf of spectra l type M5.5 (iR eid et al .||I995h 
or M6.0 (Kir kpatrick et al.lll99Ih . IFuhrmeister et al. (20051) fit- 
ted reff=2900K and \ogg - 5.5 based on a model grid with 
AFeff = 100 K and Alogg = 0.5, while Pavlenko et al. (200^ 
found Teff - 2800 ± 100 K at a fixed log ^ of 5.0 and solar metal- 
licity. From Mboi = 12.13 + 0.10 the latter determined a mass 
of 0.07-0.1 Mq and a rather young age of 0.1-0.35 Gyr, con- 
sistent the lower limit of 0.1 Gyr deduc ed from the absence of 
lithium sig natures (Maga zzu et al.lll993h. However, based on its 
kinematics. iDelfosse et al. (1998) and Mo hantv & Basril (l2003h 
classify CN Leo as intermediate young/old disk or old disk star. 

Neverth eless, CN Leo is w ell-k nown to show several signs 
of activity. iLacv et al.l (1 19761) and iGershberg & Shakhovskaial 
(1983) found high optical flare rates, howeve r, at the lowest 
flare e nergies in their samples of UV Ceti stars. iRobinson et al.l 
(Il995h analyzed UV flare events on CN Leo down to the 
microflare scale observed with the High Speed Photometer 
onboard HST. From 6 cm and 20 cm VLA observations, 
lO'Dea & McKinnonI (Il987l) could only specify a 3cr upper limit 
on its radio flux, thus CN Leo is underluminous i n this wave- 
length regime compared to other active M dwarfs. ' Giidel etal] 
( 1993) were able to detect CN Leo in the 3.5 cm and 6 cm bands 
during a flare in a simultaneous observing campaign with the 
VLA at radio wavelengths and with the ROSAT PSPC in X-rays. 
During quiescence, how ever, again only upper limits were found 
at radio wavelengths. Reale & Micelal (1199 8) discussed the X- 
ray data of the flare and derived a loop half-length of » 7 10'^ cm 
from the analysis of the flare decay. The first X-ray detec- 
tions of CN Leo ha ve akeady been reported from Einstein data 
dVaiana et al.lll98l"l : iGo lub 1983), and even during quiescence, 
it is a relatively strong X-ray source with logLx = 26.97 and 
27.01 in the ROSAT All-Sk y survey and a second po inted PSPC 
observation, respectively i lSchmitt & Liefkd |2004|) . CN Leo 
is the only star where persistent but variable optical coro- 
nal Fexm emis sion has been confirmed ( Schmitt & Wichman^ 
120011: iFuhrnieister et all 1200 4). Fe XIII is an indicator for cool 
coronal plasma (»; 1 .5 MK), while typical coronal emission mea- 
sure distributions of earhe r M dwarfs peak around 7-8 MK (e. g. 
Robrade & Schmitli l2005h . i.e., temperatures where the ioniza- 
t ion e quilibrium is shifted towards Fe xvii. Fuhrmeister et aTj 
(|2007|) found significant amounts of cooler plasma in the corona 
of CN Leo from X-ray spectra obtained with XMM-Newton, 
consistent with the Fe xiii line fluxes from simultaneous optical 
spectroscopy. In addition to its high X-ray luminosity, CN Leo 
sho ws strong Ha emission (log LHg/-Lboi - -3.89 or -3.39: 
Mo hantv & B asri 2003; IReiners & Basri 2007) and a multitud e 
of o ther chromospheric emission lines (Fuhrmeister et al.ll2004l) . 
Rei ners & Basr il (120071) measured an integrated magnetic surface 
flux of 2.4 kG, that varies by a: 100 G on timescales ranging from 
hours to days CReiners et al...2007,) . 

3. Observations and data analysis 

The giant flare on CN Leo discussed in this series of papers has 
been observed during a multiwavelength campaign on this star 
in X-rays and in the optical. This campaign consisted of six half 
nights of optical high-resolution spectroscopy monitoring with 
UVES at the VLT, and simultaneous XMM-Newton observations. 
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The first three observations have been performed on 19/20 May 
2004, and 12 and 13 December 2005; they are discussed in de- 
tail bylfuhrmeister et al. (2007). The last set of observations has 
been performed on 19/20, 21/22, and 23/24 May 2006; here we 
focus on the giant flare that occurred in the first of these three 
nights. 

XMM-Newton consists of three co-aligned X-ray telescopes, 
accompanied by the Optical Monitor OM, an optical/UV tele- 
scope that can be used with different filters in imaging or 
fast readout mode. The OM data of the flare, obtained in fast 
mode with the U ban d filter, has already been discussed by 
iFuhrmei ster et al.l 1*2008) in conjunction with the behavior of the 
optical chromospheric emission lines. XMM's X-ray telescopes 
are equipped with EPIC (European Photon Imaging Camera) 
detectors. X-ray CCDs providing medium-resolution imaging 
spectroscopy with E/AE x 20-50 and timing analysis with a 
time resolution at the subsecond level. There are two identical 
EPIC MOS detectors (MOS 1 and MOS 2), operating in the en- 
ergy range of 0.2-12.0 keV, and one EPIC PN detector, which 
covers the energy range of 0.2-15.0 keV. The PN has a higher 
sensitivity, while the MOS detectors provide better angular and 
spectral resolution; several filters and operating modes are avail- 
able for these instruments. The two X-ray telescopes with the 
MOS detectors are additionally equipped with reflective gratings 
and their corresponding CCD detectors. The Reflection Grating 
Spectrometers (RGS 1 and RGS2) provide high-resolution spec- 
troscopy in the energy range of 0.35-2.5 keV (5-38 A) with 
E/AE ranging between 200 and 800 and a spectral resolution 
of ^0.06 A FWHM which allows to resolve individual emission 
lines. All instruments are usually operated simultaneouslyQ 

In the XMM-Newton observation of 19/20 May 2006 (ObslD 
0200530501), the EPIC detectors were operated in Full Frame 
(MOS) and Large Window mode (PN) with the medium filter 
All X-ray data were reduced with the XMM-Newton Science 
Analysis System (SAS) software, version 7.0. Spectral anal- 
ysis of the EPIC data w as earned out with XSPEC V12.3 
(iDorman & Arnau dl 1200 Ih . making use of multi-temperature 
component models with coupled abundances for each compo- 
nent. The plasma models assume a collisionally-ionized, low- 
densi ty optically-thin p lasma as calculated with the APEC 
code dSmith et al.ll2001allbl) . Abundances are calculated relative 
to solar photospheric values from lAnders & Grevess3 (Il989h . 
Individual line fluxes in the RGS spectra have been mea- 
sured with the CORA program (Ness & Wichmann 2002) using 
Lorentzian line profiles with a fixed line width of 0.06 A. 

UVES (Ultra-violet and Visible Echelle Spectrograph) is 
one of the high-dispersion spectrographs of the VLT, operat- 
ing from 3000 to 1 1 000 A with a maximum spectral resolu- 
tion of llOOOCB It was used in a non-standard setup using 
dichroic beam splitter No. 2 (D1C2), with spectral coverage from 
3050 A to 3860 A in the blue arm and from 6400 A to 8190 A 
and 8400 A to 10080 A in the red arm at a spectral resolu- 
tion of =s40000. Exposure times were 1000 s in the blue arm 
and 200 s in the red arm. Here we discuss only the tempo- 
ral behavior of the optical coronal Fexiii line at 3388 A dur- 
ing the flare, which is covered by the blue arm spectra. The 
UVES spectra were reduced using the IDL-based REDUCE pack- 



CN Leo: flare light curve 



100 



' Further details on the instruments onboard XMM-Newton 
can be found in the XMM-Newton Users' Handbook, available 
at http://xmm.v ilspa.esa.eS/extemal/xmm_user_support/documentation/l 
uhb/index.html 

^ A detailed description of the UVES spectrograph is available under 
|http://www.eso.org/instruments/uves/doc/ 
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Fig. 1. EPIC PN lightcurve of large flare, corrected for pileup as 
described in the text. Time bin sizes are chosen to contain 100 
counts per bin before scaling the annular region to the total level. 
Time intervals used for spectroscopic analysis are marked. 



age (iPiskunov & Valeritill2002h . The wavelength calibration was 
carried out using Thorium- Argon spectra and resulted in an ac- 
curacy of a;0.03 A in the blue arm and »0.05 A in the red arm. 
Absolute flux calibration was carried out using the UVES master 
response curves and extinction files provided by ESO. The flux 
in the Fe xiii line has been measured with CORA. 

4. The Flare 

4.1. Timing analysis 

The XMM-Newton observation of CN Leo from 19/20 May 2006 
covers the giant flare discussed here, and two additional larger 
flares. The quasi-quiescent X-ray emission surrounding the three 
flares sums up to 20 ks, while the giant flare lasted a; 1.4 ks. The 
countrate increased from quiescent values of xO.l cts/s in the 
PN and s;0.14 cts/s in the two MOS detectors to a!65 cts/s and 
^14 cts/s, respectively, at the flare peak. Due to the long frame 
times of the EPIC imaging modes used, not only the spectral re- 
sponse but also the shape of the lightcurve is heavily affected by 
pileup. From the ratio of countrates during flare and quiescence 
from circular and annular shaped extraction regions (the latter 
excluding the innermost part of the point spread function PSF), 
we determined that the amount of pileup reached values up to 
30% for the PN at the time of maximum countrate. The effect 
on the two MOS detectors is even stronger. Figure [T] shows the 
EPIC PN lightcurve extracted from an annular extraction region 
with the innermost 150 pixels (coiTesponding to the innermost 
7.5" of the PSF) excluded. The lightcurve has been scaled by 
the ratio of the quiescent countrates in circular and annular ex- 
traction regions in order to visualize the true countrates. 

The initial rise of the few-second out burst at the beginning 
of the flare at 23:46:40 UT discussed bv ISchmitt et 31.1(120081) 
was very steep, after its short decay, the main flare continued 
to rise less steeply. However, still within less than two min- 
utes, the maximum countrate was reached. Even with the effects 
of pileup eliminated, the flare peak is flattened, resulting in a 
kind of plateau of about 3 minutes duration. The subsequent de- 
cay phase can be fitted with an exponential decay time tic of 
257±13 seconds, however, it shows slight deviations from an 
exponential shape. Especially, during the first two minutes af- 
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Fig. 2. Temporal evolution of the hardness ratio. Time bin sizes 
are chosen to contain 400 counts per bin in the total = H + S 
band. 



ter the plateau phase, the decay seems to have been faster, and 
an additional bump peaked at 0:00 UT. Less than 25 minutes 
after the initial outburst, CN Leo has returned to the quiescent 
level before the flare event. Stellar X-ray flares of that ampli- 
tude typically last much longer, i.e. from several hours up to 
days (see e. g. iKurster & Schmitl)ll996l; [Schmitt & Favatall 19991; 
iFavata et al.f l2000; Gudel et al. 2004). Integrating flie total en- 
ergy flux in the 0.2-10.0 keV band, yields a radiative loss in the 
X-ray regime of Ex ~ 4.4 -10^' erg, which can be considered 
to represent the major contribution of the total radiative loss. 
Though very energetic, the short dynamical time scales already 
suggest a very compact structure. The large flare is followed by 
a much smaller one at 0:10 UT (X-ray amplitude x2, duration 
si3 minutes), that is also visible in the optical (see Figure 1 in 
Paper I). 

In Fig. |2] we show the temporal evolution of the spectral 
hardness during the flare. We define the hardness ratio HR as 
HR - with the hard and soft bands ranging from LO- 

IS. keV and 0. 15-1 .0 keV, respectively. To avoid pileup effects, 
the annular extraction region has been used for both spectral 
bands. The shape of the HR lightcurve is fairly different from the 
normal lightcurve. Before the flare, the hardness ratio is a;-0.75, 
well consistent with what F uhrmeister et al.i (i2007l) found for the 
low-activity state of CN Leo on 19 May 2004. With the flare rise, 
the hardness ratio steeply increases, peaking at a value of s!-0.05 
about one minute before the plateau phase of the lightcurve is 
reached. At the maximum of the HR lightcurve, the contribution 
from X-ray photons above the dividing line of 1 keV therefore 
almost equalled the amount of photons below, while during qui- 
escence, the spe ctral energy distribution is by far dominated by 
lower energies. iGiidel et al.l (1200 4) found a very similar behav- 
ior for the large flare on Proxima Cen. The decay of the hardness 
ratio is interrupted by a new rise, peaking at 23:56 UT, approxi- 
mately the time when the small bump in the lightcurve starts to 
emerge. 

4.2. Plasma temperature and emission measure 

The striking changes of the spectral energy distribution are 
clearly visible in the X-ray spectra. Figure[3]compares the EPIC 
PN spectra of the quiescent state before the flare, of the flare 
rise, plateau phase and decay. The quiescent spectrum is very 
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Fig. 3. EPIC PN spectra of the quiescent state before the flare 
(black), of the rise phase (red, time intervals 1 and 2), of the 
flare peak (green, time intervals 3 and 4), and of the decay phase 
(blue, time intervals 5 to 11). 

Table 1. Fit parameters to the quiescent spectra before and after 
the flare 





before the flare 


after the flare 


Exposure time [s] 


3004 


5046 


kT, [keV] 


0.14±0.01 


0.19±0.01 


0.13±0.01 


EMi [10^" cm-''] 


16+" "3 

40+" °' 
"■^"-0.02 


0.21±0.01 


0.16+0.02 


kT2 [keV] 


0.58±0.02 


0.34±0.01 


EM2 [10'" cm-^'] 


Q 21+0.02 


0.13±0.01 


22+°-'" 


kTi [keV] 






^•^^-0.18 


EM3 [10^° cm-3] 






0.04+0.01 


red. X' 


1.05 


1.40 


1.05 


d.o.f. 


96 


149 


147 


logLx [0.2-5.0 keV] 


26.78 


26.77 


26.80 



soft, even somewhat sof ter than the spectrum from 19 May 2004. 
iFuhrmeister et al.l (12007) found CN Leo to show three differ- 
ent quasi-quiescent X-ray flux levels in May 2004 and on 1 1 
and 13 December 2006, and assumed that the two enhanced 
states in December 2006 reflect the decay phases of one or two 
long-duration flares. Our new observations now confirm the low- 
level quiescent state, which is also consistent with the previous 
ROSAT and Einstein data. 

We fitted the quiescent spectra before and after the flare 
with two-temperature component APEC models. While the tem- 
peratures and emission measures proved to be stable, rather 
large uncertainties for the individual elemental abundances oc- 
curred with the abundances treate d as free parameters. We there- 
fore chose to adopt the values of IFuhrmeister et aTl (12007 ). i. e. 
C/H = 1.49, N/H = 1.28, 0/H = 0.62, Ne/H = 1.20, Mg/H = 0.85, 
Si/H = 0.78, S/H = 1. 41, and Fe/H ^ 0.59 relativ e to the solar 
photospheric values of I Anders & Grevessd (Il989l) . The proper- 
ties of the two quiescent spectra are very similar, however, for 
the second quiescent phase the fit quality can be improved with 
a three-temperature component model. In Table[T]we summarize 
the properties of the two spectra and their associated fit parame- 
ters^ 

ISchmitt" et al.l (l2008h already showed that the flare plasma of 
the initial short-duration outburst is thermal and consistent with 
that of the following rise phase of the main flare. The large in- 



C. Liefke et al.: A giant flare on CN Leo in X-rays 



5 



crease in temperature and emission measure during the flare rise 
compared to the quiescent state is obvious from Fig.[3l the iron 
K complex clearly emerges. Enhanced line emission between 0.6 
and 1 . 1 keV sets in, when the plateau phase of the flare peak is 
reached, while during the decay both emission measure and tem- 
peratures slowly decrease. 

We divided the flare into eleven time intervals, the first six 
lasting 60 s, followed by two intervals of 120 s, and three inter- 
vals of 180 s, 240 s, and 300 s each. The first two intervals cover 
the flare rise, the next two the flattened peak, and the following 
seven intervals the different phases of the decay until the small 
flare at 0:10 UT sets in. We created spectra for each time in- 
terval and fitted them with various combinations of APEC mod- 
els, using different numbers of temperature components and sets 
of variable elemental abundances as described in the following. 
Our models always include the quiescent emission, i. e. the pa- 
rameters of the first two temperature components are fixed to the 
plasma properties of the quiescent spectrum before the flare as 
listed in Table [T] With this approach, we account for the contri- 
bution of the quiescent corona to the overall X-ray emission and 
neglect only the quiescent emission of the active region (flaring 
loop or arcade), where the flare takes place. This affects primar- 
ily the analysis of the last phases of the decay where the plasma 
properties have almost returned to the quiescent state, while the 
contribution of the quiescent emission is negligible from the flare 
rise to the middle of the decay. 

The eleven flare spectra have then been fitted with one, two, 
and three additional temperature components. The last spectrum 
however contains very few counts, considering the resulting few 
degrees of freedom, only a one-temperature component fit is ap- 
propriate. In order to compare the plasma properties of the dif- 
ferent sets of models with n additional temperature components, 
we make use of the total the emission measure EM, i. e. the sum 
of the emission measures of each temperature component / 

n+2 

EM^Y^EMi (1) 

and the flare temperature T, 
Ti ■ EMi 

r = > (2) 

^ EM 

defined as an emission measure weighted sum of the tempera- 
tures from each flare component. 

In a first approach, we fixed the elemental abundances for 
the flare plasma at the values listed above. One additional tem- 
perature component turns out to fit the last two spectra at the 
end of the decay well, while rather large residuals remain un- 
til the middle of the decay. However, already two temperature 
components give a reasonable fit for these phases of the flare. 
Even for the well-exposed rise and peak spectra, the addition of 
a third component does not provide an improvement in terms of 
X^, and models with two and three temperature components pro- 
vide consistent flare temperatures and total emission measures. 

The flare temperature peaks sharply at a!35 MK in the second 
spectrum, and immediately starts to decay. For the first six spec- 
tra, the flare-temperature lightcurve resembles the hardness-ratio 
lightcurve very well, which reflects the high-temperature contin- 
uum enhancement at higher energies. In the later phases of the 
flare, line emission becomes more important, disturbing the di- 
rect relationship between hardness ratio and temperature. The 
total emission measure shows the flat top already known from 
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Fig. 5. Temporal evolution of the iron abundance (relative to 
Anders & Grevesse (1989)). The quiescent value of 0.59±0.07 
found bv Fuhrmeister et al.l (120071) is shown for comparison. 



the lightcurve, and also the bump at 0:00 UT. Compared to the 
flare temperature, the decay is delayed by about two minutes. 

When using two additional temperature components to fit the 
flare plasma, Tt, varies between 6 and 1 1 MK, while the hot- 
ter component T4 covers a temperature range between 15 and 
45 MK, with both temperature components decreasing starting 
with the second time interval. The ratio of the two corresponding 
emission measures EMt,I EM4, however, continuously increases 
until and including the seventh time interval; the growing frac- 
tion of lower-temperature plasma shows the cooling of the flar- 
ing material. The smaller secondary peak in the hardness ra- 
tio lightcurve covered by the eighth spectrum makes EM3/EM4 
temporarily decrease again. 

4.3. Abundance variations 

Even with two and three additional temperature components, 
residuals remain between 0.8-1.1 keV and and also at ^2 keV, 
where many emission lines are located. This can be confirmed 
from the RGS spectra: Separating the RGS data into flare and 
quiescence, not only the continuum enhancement but also the 
flux increase in many emission lines in the wavelength range of 
10-17 A becomes clearly visible. Since many of the lines lo- 
cated at these wavelengths originate from highly-ionized iron, 
this may however also be a temperature effect. In order to im- 
prove the fit and to disentangle temperature and abundance vari- 
ations, we introduced single elemental abundances as additional 
free parameters to the sequence of fits. Particularly, we tested 
iron, neon, silicon and oxygen. 

With the iron abundance as a free parameter, the maximum 
flare temperature is reduced to 32 MK; a strong enhancement in 
iron during the flare compensates part of the flux of the iron K 
complex, which serves as a major temperature indicator The to- 
tal emission measure only slightly decreases. However, the gen- 
eral behavior of the one-, two-, and three-temperature compo- 
nent fits remains unchanged, so the shape of both the flare tem- 
perature and the total emission measure lightcurves is preserved, 
as shown in Fig |4] With the iron abundance as a free parameter, 
the fit quality clearly improves. Table |2] lists the fit parameters 
for the optimum combination of two flare components for the 
first ten spectra and one component for the last spectrum. 
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Fig. 4. Temporal evolution of flare temperature (left) and total emission measure (right) of the flare plasma obtained from a two- 
temperature component fit with variable Fe abundance. 



Table 2. Model parameters for the flare spectra with variable Fe abundance 



Spect. 


Exp. 


kT, 


EM3 kT4 EM4 


Fe logLx "^^^-X^ d.o.f. 


T EM 


No. 


time [s] 


[keV] 


[10™ cm-3] [keV] [10™ cm-^] 


[0.2-10.0 keV] 


[MK] [10™ cm-3] 



-1 



l-OS^pf, 28.90 0.99 91 27.97+|'» 43.58+ 

■98 1 00+0:21 TQ m 1 1 m I"! in+lOB 



2 60 0.84+1;™ 6.99+!-| 3.06+"-J^ 52.11!||° 1-38!o;q 29.07 1.22 130 32.50+||^° 59.10!?-| 

3 60 0.75+0.02 9.57+if2 2.59+|T^ 50.60!j-| l-SS^j^ 29.08 1.23 141 26.69+|| 60.18!^]^ 

4 60 0.73±0.03 10.97!|-^^ 2.34+||l5 48.06!^-°* l-O?^!! 29.05 1.09 132 23.66+if 59.03!|°^ 

5 60 0.75+i!!]2 14.46+:]^ 2.08+!!" 37.13+^:* 0.89+™ 28.99 1.10 121 19.84+[" 51.59+^-2' 

6 60 0.71+11™ 10.67+;:'^ 2.02+||-;| 28.15+^"*? 0.72+|j-j* 28.84 1.02 91 19.27+i-[» 38.82+^™ 

7 120 0.71+'j<]3 6.32+^]^ l-91!o20 H-^S^^qJ 0.75+;J;3 28.52 1.26 91 17.27+|| 18.07+f^» 

8 120 0.61+1111^ 1.64+»39 1.53+o[^ 7.05+!!| 0.85+!![^ 28.20 0.89 46 15.75+|| 8.69+»|^ 

9 180 0.58!!J°^ l-64+!!t? l-31!o!? 3-33!™ 0.59!||;* 27.94 0.91 40 12.42+[{^ 4.98!J!-?? 

10 240 0.55+!;-?^ 1-24!|;-^] l-35!oi8 l-13!o32 0-52!on 27.62 0.95 28 10.85+|-'^ 2.37!;;-^;; 

11 300 0.71+°')^ 0.36±0.11 0.49+°-^^ 26.84 0.80 10 8.23;|-/,^ 0.36±0.11 



Figure |5] shows the corresponding time serie s of the 
iron abundance relative to Anders & Grevesse' (119891) with the 
Icr confidence range o f the quiescent abundance found by 
iFuhrmeister et al.l (I2007D marked as a shaded area. The Fe abun- 
dance increases from a clearly subsolar level of 0.59±0.07 dur- 
ing quiescence by more than a factor of two to a maximum 
value of 1.38±0.20, when the flare peak is reached. It then de- 
cays even faster than the flare temperature; consistency with the 
quiescent value is reached again with the 6th spectrum. A re- 
increase at 23:55 UT in the eighth spectrum may not be signifi- 
cant, however, it coincides with the second rise of the hardness 
ratio lightcurve. By the end of the decay, the iron abundance has 
eventually dropped to the quies cent value again . A ver y simi- 
lar behavior has been found by iFavata & Schiiiiti (1 19991) for the 
overall metallicity during a long-duration giant flare on Algol 
observed with BeppoSAX. 

Even with the Fe abundance as a free parameter, the second 
and the third spectrum cannot be fitted with a statistically accept- 
able value xled- Residuals remain mainly around 1 keV in the 
second and between 0.7 and 1.2 keV as well as at 2 kev in the 
third time interval, respectively. There are quite a few emission 
lines at these energies, not only from iron but also from oxygen, 
neon and silicon. This indicates that further abundances may de- 
viate from the quiescent values. When the neon abundance is set 
free, it seems to be generally enhanced compared to the value 



of 1.20+0.16 from IFuhrmeister et all (l2007 l). the Ne abundance 
does however not give a clear pattern like the iron abundance in 
Fig. |5]but shows a large scatter, and the uncertainties are rather 
large. Most values are on the Icr level consistent with the old qui- 
escent value. A free neon abundance therefore only marginally 
improves the fit, most of the residuals remain, with flare temper- 
ature and emission measure consistent stable at the initial val- 
ues. Silicon shows a similar behavior, with its abundance values 
scattering around the quiescent value of 0.78+0.33. Again, tem- 
perature and emission measure stay at the values obtained from 
the fitting sequence with all abundances fixed. The same proce- 
dure with oxygen tends to result in unreasonable fits with either 
the oxygen abundance or at least one of the temperature com- 
ponents diverging or approaching zero. We also combined two 
or more free elemental abundances, i. e. Fe and Ne, Fe and Si, 
as well as Fe, Ne and Si, but none of these combinations pro- 
vided further improvement in the quality of the fit compared to 
Fe set free only. These fits even start to diverge, especially for 
the less well-exposed last two spectra with few degrees of free- 
dom. Additionally, the residuals in the fits of the second and the 
third spectrum cannot be reduced significantly. 

Iron therefore seems to be the only element, where a clear 
variation of the elemental abundance during the flare and com- 
pared to the quiescent state can be observed. In the common pic- 
ture of stellar flares, the evidence of changes in the iron abun- 
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dance indicates that the flare plasma, i. e. fresh material evap- 
orated from photosphere and chromosphere, shows a different 
composition with a higher iron abundance compared to the or- 
dinary coronal plasma. Since CN Leo's quiesc ent corona shows 
the invers e FIP effect (iFuhrmeister et al.l2007h . this confirms the 
finding of lNordon & Behad (l2008h that IFIP-biased coronae tend 
to show a FIP bias during flares compared to their quiescent 
state. However, the Fe abundance seems to decrease much faster 
than temperature and emission measure of the flare plasma, and 
the iron abundance has returned to values consistent with the 
quiescent level when the X-ray emission of CN Leo is still by 
far dominated by the flare. The time-resolved evolution of the 
iron abundance in the flare plasma suggests either the existence 
of some kind of iron-depleting mechanism in the flare plasma, 
which would be rather difficult to explain, or the possibility that 
different parts of the lower atmosphere with different abundance 
levels were affected during the early and later stages of the evap- 
oration process. 



Finally we would like to address the possibility that system- 
atic errors in our fit models actually mimic an enhanced iron 
abundance. At the very hot plasma temperatures as encountered 
during a flare, the iron abundance fit is affected by the flux of 
iron lines from lower ionization stages around 0.8-1.2 keV and 
- more importantly - by the flux of the emission complex at 
6.7 keV, mainly originating from iron in ionization stages from 
Fe XXIII to Fe xxv. An incorrect temperature model would lead to 
an incorrect prediction of the emission both in the 0.8-1.2 keV 
energy range and in the 6.7 keV line complex, which would 
then have to be compensated by changing the iron abundance. 
Whether this is the case or not is hard to tell, yet a detailed study 
undertaken by dFavata & Schmittll 1999 1) suggests that this is un- 
likely to happen. Further, Fe Ka fluorescence at 6.4 keV, excited 
in the neutral photospheric material by intense high-energy ra- 
diation >7.1 keV from the flare plasma, may also contribute. Fe 
Ka fluorescence has been observe d during strong stellar flares 
(e.g. lOsten etani2007HTesta et al.'2008). The Ka excitation de- 
pends on the photospheric iron abundance and on the inclination 
angle and he ight of the X-ray emitting flare plasma above the 



CN Leo OVIl triplet quiescence 



photo sphere dTesta et alJ2008l : lErcolano etal .l2008HDrake et alJ 
|2008|) . A compact loop structure as probably responsible for the 
short decay time of the flare on CN Leo is expected to produce 
rather strong fluorescence, yet no excess emission at 6.4 keV 
is visible, especially during the second and third time intervals, 
where the spectral fits yield the largest Fe abundance values. All 
eleven spectra are well consistent with models without any Fe 
Ka fluorescence emission at 6.4 keV. A possible explanation for 
this finding is an unfavorable viewing angle. At any rate, the ex- 
pected Fe Ka line features are relatively weak, they can be spec- 
trally separated from the 6.7 keV complex and should therefore 
not lead to incorrect Fe abundance estimates. Deviations from 
ionizatio n equilibrium may also affect the strength of the 6.7 keV 
complex: iReale & Orlandol ([200 8) have shown that the effect is 
significant for flares with heat pulse durations less than a minute. 
However, it is difficult to quantify this effect for the CN Leo flare, 
also such deviations from ionization equilibrium are strongest 
at flare onset, while ionization equilibrium is established rather 
quickly with time in a high density environment as encountered 
with the flaring plasma in CN Leo, so non-equilibrium effects 
appear to be an unlikely cause for the observed changes in iron 
abundance. Thus in summary, no clear systematic effects are 
known that would mimic an enhanced iron abundance during 
the flare evolution. 




21.5 22.0 
wovelength [ A ] 



CN Leo OVIl triplet flare 




21.0 21.5 22.0 22.5 

wovelengtti [ A ] 

Fig. 6. The O vii triplet in the RGS spectra of CN Leo during qui- 
escent (top) and flaring (bottom) state. The flare spectrum covers 
the same time interval as the eleven EPIC spectra, the quiescent 
data include the time intervals before and after the giant flare, 
with two further flares excluded. 



4.4. Densities 

The RGS spectra allow to investigate the electron densities of 
the coronal plasma from the density-sensitive ratio of the for- 
bidden and intercombination lines of helium-like triplets. With 
this method, Giidel et al.l ( 120021) found the density of the coro- 
nal plasma as traced by the Ovn and Neix triplets to vary by 
more than an order of magnitude during a giant flare on Proxima 
Centauri, which is the only unambiguous measurement of in- 
creased coronal densities during a stellar flare. While the two 
flares were comparable in maximum countrate, the duration of 
the CN Leo flare was much shorter, so that, in contrast to the 
Proxima Centauri flare and to the CN Leo EPIC data, the CN Leo 
RGS data can only be separated into quiescent state and flare; the 
RGS signal is too low for further sub-division. 

In Fig. |6] we show the RGS 1 spectrum centered on the O vii 
triplet during the 20 ks of accumulated quiescence and during 
the flare together with the best fits to the three triplet lines r (res- 
onance), / (intercombination), and / (forbidden) provided by the 
CORA program. The measured line counts, deduced // / ratios 
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Table 3. Measured line counts, // / ratios and deduced coronal 
densities from the O vn triplet in the RGS spectra of quiescence 
and giant flare in Fig.|6] 



Line 


Wavelength [A] 


Quiescence 


Flare 


vn r 


21.60 


30.2±6.8 


25.4±5.9 


vn ( 


21.81 


9.4±4.5 


14.0±4.7 


Ovn/ 


22.10 


27.1±6.3 


1.2+2.0 


f/i 




2.88±1.53 


09+" ' 


log lie 




< 10.77 


>11.70 



and electron densities are listed in Table [3] In order to convert 
the measured /// ratios to densities, we used the relation 



with the critical density A^^. and the low-density limit Rq, 
where we adopt e d valu es of 3.1 ■ 10"^ cm"^ and 3.95 from 
iPradhan & Shulll (Il98lh . During quiescence, the /// ratio is 
well consistent with the low-density limit, which is in good 
agreement with the values found for CN Leo in May 2004 
by Fuhrmeister et al. (2007). In the flare spectrum, the forbid- 
den line has almost vanished, while the intercombination line 
is stronger The electron density therefore clearly deviates from 
the quiescent value of rig - 1.2^j 2 ■ 10'" cm"-' during the 
flare, and is definitely incompatible with the low-density limit. 
The average density during the flare has a nominal value of 
He = 1-4 ■ 10'^ cm"^. Taking into account the Icr error on the 
line fluxes, a lower limit of > 4.9 ■ 10" cm""*, i.e. at least 
an order of magnitude higher than the quiescent value, can be 
derived. 

/// ratio and electron density obtained from the RGS flare 
spectrum represent mean values for the whole flare. In order 
to give a rough estimate of the peak flare density, one can de- 
rive a modifying scaling factor of xl.2 based on the square 
root of the maximum and mean total emission measure values 
from the EPIC data and assuming a constant plasma volume 
during the flare. The flux in the forbidden and intercombination 
lines measured during the whole flare is therefore dominated by 
the high countrates at the flare peak, and with a lower limit of 
rie > 5.9 ■ 10" cm"^ the maximum density is only slightly larger 
than the average value. 

With a peak formation temperature of ^2 MK, the Ovii 
triplet traces only the cooler component of the coronal and flare 
plasma, and the density sensitivity of its /// ratio is optimal be- 
tween 10'" and 10" cm"-*. As the O vii triplet suggests densities 
probably even higher than 10'^ cm"-', it would be desirable to 
probe higher densities and the higher-temperature plasma with 
the helium-like triplets of heavier ions. Unfortunately, the flux 
in the RGS2 spectrum in the Neix triplet, where the density- 
sensitive range has shifted to 10 "-10'^ cm"-' and the peak for- 
mation temperature has increased to ^4 MK, is too low and its 
contamination with highly-ionized iron during the flare is too 
strong to give any conclusive result. Mg xi and Si xiii cover even 
higher densities up to lO'"* cm"^ and temperatures up to 10 MK, 
but in the two RGS spectra, the resonance lines of both triplets 
are only barely visible. A density survey of the flare at different 
plasma temperatures is therefore not possible. 



4.5. The optical coronal Fexiii line 

With the RGS, the wavelength range a!6-37 A can be observed; 
including individual emission lines with peak formation tem- 
peratures from a: 1.5-2 MK in the low-temperature range (Cvi, 
Ovii) up to 15 MK (Sixiv). The optical coronal Fexiii line at 
3388 A included in the blue-arm UVES spectra provides addi- 
tional time-resolved information on the lowest coronal tempera- 
tures (peak formation temperature ^1.5 MK). 

Compared to the previous observations in December 2005, 
the Fe xiii line is rather weak in all three observations in May 
2006. The measured line fluxes have large uncertainties, some- 
times the line is not even clearly detected. The measured FWHM 
is also smaller than previously observed. Unfortunately, no 
Fe XIII line data is available for comparison from the observation 
in May 2004, where CN Leo's quiescent X-ray emission was at a 
lower level, similar to the quiescent phase outside the large flare. 
However , it is reasonable to assume that the overall reduction in 
emission measure, which affects the lower temperatures where 
Fe XIII is formed as well, is responsible for the lower line fluxes. 

The most significant detections are obtained after the giant 
flare and after a second larger X-ray flare. A third larger X-ray 
flare is not covered by the optical observations anymore. The 
Fexiii line is not visible in the pre-flare spectrum and in the 
two flare spectra; but during the flare, a chromospheric Co i line 
emerges at the same wavelength position. 

Even during quiescence, the bulk of iron in the coronal 
plasma should exist in the form of neon-like Fe xvii, with only 
minor contributions of the cooler Fe xiii. With the onset of the 
flare, the temperature of the X-ray emitting plasma shifts to even 
higher values, and therefore the ionization balance should be 
dominated by far by higher ionization stages. Even with an in- 
creased iron abundance, Fe xiii can be expected to be depleted 
until the end of the decay, when the plasma has cooled down 
again. The excess emission measure then allows the detection of 
the Fexiii, which is consistent with the observations. 



5. Loop modeling 

Direct imaging demonstrates the great complexity of flaring ac- 
tive regions on the Sun. In order to obtain a consistent model 
describing the physical properties and evolution of the involved 
structures, a stellar flare can be reduced to a sequence of pro- 
cesses taking place in a simplified geometry. Such a gener- 
alized description starts with an initial heat pulse in a single 
coronal loop and efficient conduction until the whole loop has 
reached the maximum temperature. Lower atmospheric layers 
are strongly heated and chromospheric evaporation fills the loop 
with denser material. Conductive cooling sets in while the filling 
of the loop may still go on. Radiative cooling dominates as soon 
as the maximum density is reached, and loop depletion is initi- 
ated. Cooling and loop depletion may be decelerated by potential 
residual heating. The transitions between these processes can be 
described in terms of equilibrium conditions, but their possible 
overlap and interplay complicate the modeling. However, over 
the last decades, model approaches for flares have evolved from 
scaling laws valid only for (quasi)static conditions or simple an- 
alytical estimates to complex hydrodynamic loop models taking 
into account various initial conditions. These models have been 
successfully applied to a wide variety of solar and stellar flares. 
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5.1. The decay phase 

A basic finding for solar and stellar flares is that the decay time 
of a flare scales with the length of the flaring loop. In order to 
quantify this effect, Serio et al. (1991) developed an analytic ap- 
proximation relating the thermodynamic decay time t,/,, the loop 
half length L and the loop top maximum temperature To of a 
flaring loop, based on the underlying set of hydrodynamic equa- 
tions and assuming semicircular loops with constant cross sec- 
tion cooling from equilibrium conditions and uniformly heated 
by an initial heat pulse without any furt her heating during the 
decay. This approach has been refined by i Reale et aP (Il993h to 
account for loops comparable or larger than the pressure scale 
height h at the loop top 



h = 



fig 



(4) 



where kg is the Boltzmann constant, fi is the effective mass per 
particle and g the surface gravity of the star. The pressure scale 
height for the giant flare is great e r than 100 000 km. 

In their models, iReale et aP (Il993i) also consider the effect 
of a gradually decaying heating functio n, and finall y , for l oops 
smaller than the pressure scale height, IReale et alJ (Il997h de- 
rived and tested an empirical expression to determine the flar- 
ing loop length including the effect of sustained heating during 
the decay that uses the slope ( of the flare decay in the density- 
temperature plane to evaluate the amount of sustained heating 
during the decay: 



tlc 



^/T'o 



aF(0 



or 



120 F(0 



(5) 



lO""* cm-i s K'/2. 



where a = 3.7 

The term F(0 describes the ratio of observed and thermo- 
dynamic decay time TidTth as a function of the ratio of tem- 
perature and density decay ^. The higher F{^), the larger is the 
amount of prolonged heating during the decay, while F{^) « 1 
would indicate the absence of additional heating. F{^) is mono- 
tonically decreasing and can be approximated with analytical 
functions; here we use the hyperbolic form 



Tth (,-(,a 



+ qa 



(6) 



The coefficients q,, ^a, and q^, depend on the energy response 
of the instrument used. For XMM/EPIC, the v alues a r e Cg - 
0.51, = 0.35, and = 1.36 according to iReald (l2007h . 
Additionally, the observed maximum of the flare temperature 
ro.obs, which is an averaged value over the whole loop, must be 
corrected with 



7^0=^ Tl^ 



To,i^^-m-' -Tl 



(7) 



in order to obtain the true Tq independently of the detector 
The values o f the coeffic ients are f - 0.13 and rj - 1.16 for 
XMM/EPIC (iRealell2007h . With the maximum flare temperature 
found in the second time interval, we determined Tq -j- 5.6+0.4. 

Figure|7]shows the evolution of the giant flare in the density- 
temperature plane, where we use VeM as a proxy for the den- 
sity. We fitted the slope f of the decay, i. e. with the data from 
first three spectra excluded, to ^ = 0.40 + 0.05. This value is 
rather small, and on a Icr level consistent with the minimum 
value of 0.35, which means that the limit of validity for this 
method is reached. For very small the corrective term F{() 



7.6 



7.4 



CP 7.2 



7.0 



6.8 



= 257 s 



( = 0.40 




L = 490013900 km 



24.8 25.0 25.2 25.4 25.6 25.8 26.0 
log VEM 

Fig. 7. Flare evolution in the density-temperature plane. 



that accounts for the presence of heating during the decay be- 
comes very high and also rather uncertain, as the hyperbola ap- 
proaches infinity. For ^ - 0.40, we obtain a formal value of 
F{^) - ticItiIj - 10.7; the thermodynamic decay time would 
therefore lie in the range of ^ 25 seconds, which is extremely 
short and hints to very effective cooling. The size of the involved 
loop structure is thus very small, the resulting loop half length 
is L = 4900 ± 3900 km. The uncertainty in L is dominated by 
far by the uncertainty in ^ in this flat part of the hyperbola of 
F(^). Effectively, we can therefore only give an upper limit for 
L, which can be set to ^ 9000 km using the Icr error 

The original approach of i Reale e t al.l ( |l997 l) is based on hy- 
drostatic and energy equilibrium at the beginning of the decay. 
This would result in a simultaneous onset of the decay of temper- 
ature and density, which is apparently not the case for the giant 
flare on CN Leo. In order to account for this, the flare tempera- 
ture at the time of maxim um de nsity Tu should be used instead 
of Tq according to ReajJ (l2007h . otherwise the loop length will 
be overestimated. However, considering that the dependency of 
L on the square root of Tq is comparably weak, we keep at the 
original value for the upper limit. Additionally, the exact point 
in time of maximum density is difficult to assess from flat top of 
the emission measure lightcurve, while the temperature already 
decreases significantly. 



5.2. Flare rise and peak 



Following the approach of lReald ( l2007h . the duration of the flare 
rise and a possible delay between the decay of temperature and 
density can also be used to characterize the dimensions of a flar- 
ing loop. Both time scales can be traced back to equilibrium 
states of heat input and cooling. The occurrence of such a de- 
lay, as it is also observed for the giant flare on CN Leo, indicates 
an initial heat pulse too short to bring the flare plasma to equi- 
librium conditions, again pointing toward short time scales and 
a compact structure. 

The maximum density is reached when conductive and ra- 
diative cooling equal. This gives 



(8) 



Lg -3iff -J To J fM. 
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as an estimate for the loop length from time tu at which the 
maximum density occurs and the ratio ij/ of the maximum tem- 
perature T{) and the temperature at maximum density Tm'- 



1 M 



(9) 



From the time AJo-m - tu — to between the start of the tempera- 
ture decay and the density maximum, i. e. the start of the density 
decay, one obtains 



adopting the conductive cooling time. 

These two independent approaches are particularly useful, if 
only the earlier phases of the flare are observed, i. e. Eq.|8]can be 
applied if the complete rise phase of the flare is available until 
the density maximum is reached, while bigger parts of the decay 
are not required. Eq. [10] can be used if additionally the quies- 
cent phase before the flare and the flare onset itself are missing. 
When the flare is completely covered by the observations, the 
two methods yield independent estimates of the loop half length 
which can be cross-checked. 

Both Eqns. [8] and [TOl depend strongly on if/ and also on the 
two onset times of the temperature and density decay, which 
have to be determined as accurately as possible. Using the time 
intervals from Sect. 14.21 the temperature maximum occurs dur- 
ing the second spectrum, which would result in fo.spec = 90+30 s, 
i. e. a rather large uncertainty due to the length of the time inter- 
val compared with the rise time itself. It is, however, not rea- 
sonable to divide the sequence any further. The hardness ratio 
lightcurve from Fig. |2] provides a better time resolution, assum- 
ing that especially during the early phases of the flare, the hard- 
ness ratio represents a good temperature indicator. The hardness 
ratio lightcurve undergoes major changes during the first time in- 
terval; it peaks already at fo.HR - 47+5 s, earlier than the spectral 
fit sequence, that provides only the average temperature of each 
interval. We would therefore prefer to use the hardness ratio as 
a temperature tracer, but subsequently, the true maximum tem- 
perature Tq will be larger than the 56 MK derived from the flare 
temperature of the second time interval, and it is difficult to give 
an accurate estimate. The difference between the two possibili- 
ties to choose for fo will at any rate strongly affect the results of 
the loop length estimate. 

It is even more difficult to define the maximum of the den- 
sity. The emission measure, which we have used as the density 
indicator in Sect. 15.11 peaks in the third time interval. However, 
due to the flattened shape of the emission measure lightcurve, the 
uncertainty in fM is probably even larger than the corresponding 
value of 150 ± 30 s. Consistent results are obtained when us- 
ing a smoothed version of the countrate lightcurve, which would 
nevert heless result in an underestimate of according to lReald 
(I2007h . Tm is therefore determined from the flare temperature in 
the third time interval. 

Using the numbers obtained from the analysis of the time in- 
tervals only, we obtain L = 25 000 ± 18 400 km from Eq.fTOl 
When sticking to To from the second time interval but in com- 
bination with fo from the hardness ratio lightcurve, L increases 
to 43 000 + 15 300 km, and a higher value for To would even 
result in a further increase, with the corresponding value being 
unknown. Eq.[8]gives L = 16 700 + 4900 km, closer to the value 
from the upper limit from the analysis of the decay. Considering 
the large error bars, the results from Eqns.lSlfTOl and[8]are con- 
sistent at least on the 2cr level. The overall uncertainties are dom- 
inated by the uncertainties of the two time scales fo and fM: For 



the results obtained from the flare rise from Eq. [8] the uncer- 
tainty in fM makes up 53% of the uncertainty in L, while To and 
Tm share 33% and 13% respectively. When the time intervals are 
used to determine fo in Eq.[TOl Afo-M accounts even for 91% of 
the uncertainty in L. Using fo from the hardness ratio lightcurve, 
this reduces to 63%. As Eq.[TO]yields indeterminate results, Eq.[8] 
should be preferred. 

6. Discussion 



(10) 6.1. implications of the modeling 



The giant flare on CN Leo has a very short exponential decay 
time of only a few minutes compared to other flares with sim- 
ilar amplitudes in X-rays, which typically last for hours. This 
holds even if one considers only flare s on active M dwarfs (e. g. 
iFavata et al.] l2000': Gudel et al.l l2002h . Such flares mostly origi- 
nate from large loop structures with dimensions of at least a sig- 
nificant fraction of the size of the star itself. The short duration of 
the CN Leo flare on the other hand points at short cooling times 
and therefore a rather compact involved structure. Respective es- 
timates assuming a single flaring loop have confirmed this pic- 
ture, the resulting loop half length is of the order of about a tenth 
of the stellar radius or even smaller. 

We can assess the emitting plasma volume V via 



EM 



(11) 



and obtain the area A of the loop footpoint on the steUar surface 



(12) 



approximating the loop with an unbent tube with a circular base 
area and constant cross-section over the loop length. In order to 
compute V and A, we can use the electron density derived from 
the helium-like O vii triplet of the flare spectrum in Sect. 14.41 
and the total emission measure averaged over the eleven time 
intervals from Sect. |42] EM - 1.8 ■ 10^' cm ^ and obtain 
V < 7.4 ■ 10^^ cm^ from the lower limit of the density, or 
y = 9.1 ■ 10^'' cm-' using the nominal value. This results in values 
of 1.3-80 ■ 10'^ cm^ for A, depending on the value used for L. 
Larger loop lengths, as they are provided by the data from the 
flare rise and peak, give the smallest footpoint area. A good ref- 
erence point is A « 1 ■ 10'** cm^, which would correspond to a 
short loop with a length close to the nominal value derived from 
the decay (i. e. s!4500 km) at the nominal electron density value 
of «£ = 1 -4 • 10'2 cm"^, or a longer loop with L a;35 000 km at a 
lower density of 5 • lO" cm"^. 

Another possibility to estimate independently from the 
emission measure is to make use of the scahng laws of 
iRosner et al ] (ll978h 



(13) 



To = 1400 Vpl 



which assume equilibrium conditions and the plasma equation 
of state p - IriekBT: 



no = 1-3- 10*"-^ 



il. 

L 



j2 

no.io = 13-;— 



(14) 



This results in density values between - 3.9 ■ lO" cm"-' and 
«£ = 2.8 ■ 10'2 cm"-' for the largest and smallest possible loop 
lengths respectively, which matches the maximum density esti- 
mate from Sect. l4.4l quite well. The maximum density «m should 
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be somewhat higher than no- However, due to the large uncer- 
tainties in L, the derived densities cannot be fixed to an order of 
magnitude. Similarly, footpoint areas between A = 7.7 ■ 10'^ cm-' 
and A = 3.2 ■ 10'** cm-' are obtained, but again with large uncer- 
tainties on individual values. In addition, it is unclear to which 
level this particular flare follows the established scaling laws at 
all, as the initial assumptions like equilibrium or constant flaring 
volume are not necessarily fulfilled. 

In order to model the short initial burst, Schmitt et alj (l2008h 
assumed a single loop with a half length of 20 000 km. Along the 
loop, they obtain densities of 1-10 ■ 10" cm"^ after it is filled 
with evaporated material. The temperature in the loop model 
reaches maximum values around 4- 10^ K, both numbers are well 
co nsistent with what we obta in also for the major flare event. 

lAschwanden et al] (12008 ) compiled statistical correlations 
between peak temperature Tq, peak emission measure EM, X- 
ray luminosity Lx, total X-ray-radiated energy Ex, and flare du- 
ration Tf for a large sample of observed stellar flares. They find 
the following relations as a function of the maximum tempera- 
ture Tq, 

EMiTo) = 10^"-^ ■ T^-j^^-"^ (15) 

T/(ro) = lo^-^ ■ r^f (16) 

Lx(To) = lO^'-' ■ T*-^^°-'' (17) 

ExiTo) = 10^"-^ ■ Tg^j^"-^ (18) 

however, the correlation between Tf and (ro) in Eq. [16] is only 
marginally significant due to large scatter in the data. Our mea- 
surements from giant flare on CN Leo clearly deviate from the 
relation between peak temperature and emission measure, the 
observed maximum emission measure is more than an order of 
magnitude lower than the value provided by Eq. [15] While the 
densities we observe seem to be common for such a huge flare, 
the flaring structure is rather compact, i. e. the flaring volume is 
very small, which apparently causes the low emission measure. 
Eq. [16] predicts a flare duration of about 45 minutes, more than 
two times longer than observed. However, as the correlation be- 
tween Tf and Tq is not that strong, this matches the data regime 
from the flare sample. As a consequence of the deviation of the 
emission measure, Eq.[T7]overestimates the X-ray luminosity at 
the flare peak also by an order of magnitude. The total radiative 
loss in X-rays from Eq. [18] is too high even by two orders of 
magnitude as a result of both the overestimated emission mea- 
sure and flare duration. 

6.2. Single loop or arcade? 

The Sun shows us that the model of an accumulation of iso- 
lated single loops oversimplifies a stellar corona. Instead, loops 
of different sizes often group to complex structures of arcades, 
especially in the active regions where flares originate. Such 
loop structures are non-static; individual loops resize, their foot- 
points move, and they realign in magnetic reconnection events. 
Unfortunately, we cannot even resolve individual active regions 
directly on stars other than the Sun. Doppler imaging is a com- 
mon indirect method to reconstruct photospheric spots or promi- 
nences in the chro mosphere of fast ro tators spectroscopically in 
the optical (see e. g IStrassmeieil l2009) . Lightcurve-based eclipse 
mapping techniques pr ovide an addit i onal ap proach, as success- 
fully demonstrated by IWolter et al.l ( l2009l) . However, the so- 
lutions of such image reconstruction methods are usually not 



unique, becoming more and more ambiguous the better the tar- 
geted spatial resolution is. Additionally, such observations are 
very rare in the X-ra y regime, cf. Giidel et al. (2001, 2003), and 
ISchmittet al.l(l2003h . 

There are however other clear indications for the substructur- 
ing of active regions in stellar coronae. Flare lightcurves often 
deviate from the simple scheme of rise and exponential decay, 
and instead show multiple peaks or v ariations in the decay time, 
like a double exponential decay (e. g. lWargelin et al ] l2008l) . The 
lightcurve of the giant flare on CN Leo shows several irregu- 
larities like the initial outburst and the flat peak. Its decay de- 
viates from an exponential shape, and there is the small bump 
in the lightcurve which is apparently related to a second peak 
in the hardness ratio lightcurve. This event obviously interrupts 
the otherwise continuous cooling of the flare plasma. Such fea- 
tures suggest that multiple components of the active region from 
which the whole flare event originates are successively involved, 
each of them with its own properties like shape and size, tem- 
perature and density, and hence cooling time. 

Further evidence for flare substructure can be obtained from 
the spectral analysis of the flaring plasma. A strong correction 
F((), i.e. when heating dominates the decay of the flare, points 
to an arcade of loops involved, as observed for t wo-ribb on flares 
on the Sun; this was first pointed out by Realei (l2002h . For the 
giant flare, F{0 is very high. We therefore interpret the whole 
event in terms of a flaring arcade. In this scenario, the flat top 
of the lightcurve can be explained by the superposition of sev- 
eral loops flaring consecutively. Later during the global decay, 
further loops not affected so far became involved and caused the 
small bump. This flare series was triggered by the short initial 
outb urst, which could be localized in an almost isolated single 
loop dSchmitt et al.ll2008h . If we assume that the flat peak of the 
main flare is composed of a sequence of unresolved individual 
peaks with similar luminosity, the true values of and in partic- 
ular Afo-M would be lower while t^c will stay roughly the same, 
which would in turn reduce L from Eqns. [8] and [TO] i.e. bring 
the results from the analysis of rise and peak closer to the upper 
limit obtained from the decay. 

A theoretical model framework for solar two-ribbon flares 
has been developed by iKopp & Polettd ([1984) and - with re- 
stri ctions - s uccessfully applied also to stellar flares, see e.g. 
Pol etto et al.l (Il988) for the modeling of two long-duration flares 
on the active M dwarfs EQ Peg and Proxima Centauri observed 
with Einste i n and EXOSAT. The single-loop model approach of 
Re a le et a D (fT997l) can be extended to multi-loop structures if 
different phases of the flare are analyzed separately. iReale et aP 
(2004) investigated the evolution of the long-duration giant flare 
on Proxima Cen observed with XMM-Newton. In their best-fit 
model, the X-ray emission is dominated by a single loop for the 
first part of the flare, while later on during the decay an addi- 
tional arcade of loops with approximately the same length as the 
initial single loop emerges. This picture is supported by the anal- 
ysis of spat iaUy resolved observations of strong solar two-ribbon 
flares (e.g. lAschwanden & Alexandeiil2001,) . For the giant flare 
on CN Leo, the primary single loop was causing the 2-second 
initial outburst. Due to its short duration it is impossible to an- 
alyze this very first part of the flare in more detail in order to 
obtain flare temperature and loop length. We can, however, con- 
clude that even if the whole flare event is treated as originating 
from an arcade-like structure, it would be very compact, prob- 
ably even smaller in terms of loop length than a single loop. 
Reale et al. (2004) found a loop length larger by 30% if the flare 
on Proxima Cen is interpreted with Eqns.|5]-[7]instead of with de- 
tailed hydrodynamic loop modeling. For a single loop, the flar- 
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ing volume scales with L, while an arcade of loops implicates 
y oc L^. For solar flares, an empirical correlation V oc is ob- 
tained (jAschwanden et al. 2008), reflecting the observation that 
arcades are not uniformly filled but consist of individual loops. 

6.3. The multiwavelength context 

Plasma temperature, density, and optical thickness of the stel- 
lar atmosphere change by orders of magnitude, when going 
from the photosphere via chromosphere and transition region 
up into the corona. According to its spectral type, i.e. the ef- 
fective temperature, the major contribution to the photospheric 
spectrum of an M dwarf like CN Leo is emitted the infrared, 
with a much smaller proportion observed in the optical. The 
chromosphere adds line emission, mainly in the UV and in the 
optical. Line formation at these wavelengths is a complex pro- 
cess. Photospheric absorption lines may turn into emission in 
the chromosphere, and the strength and shape of chromospheric 
lines is extremely sensitive to the atmospheric conditions. A so- 
phisticated modeling of M dwarf chromospheres is challenging, 
and current mode ls are still far from providing a complete de- 
scription (see e.g. iFuhrmeister et al.|[2005i) . Nevertheless, as for 
each line the formation takes place under specific physical con- 
ditions which basically translate into certain heights above the 
the "surface", diff'erent lines probe the structure of the stellar at- 
mosphere. Transition region and corona emit mainly in X-rays 
and in the EUV. Complemented by non-thermal processes like 
e.g. radio gyrosynchrotron emission, the different atmospheric 
layers can therefore be studied over a broad wavelength range. 

Solar and stellar flares can affect all parts of the atmosphere, 
with the strongest impact observed on chromosphere and corona. 
However, as on M dwarfs the quiescent photospheric continuum 
is low in the optical, photospheric flare continuum emission, 
that accompanies strong events like the giant flare on CN Leo 
discussed in our series of papers, is much more obvious com- 
pared to the Sun, where white light flares are rarely observed. 
When integrating over the full disk, they would go unnoticed in 
optical broad band measurements. In Paper I, we analyzed the 
observed optical continuum emission during the CN Leo flare. 
We found that the photospheric material was rapidly heated to 
temperatures probably exceeding 20000 K at the onset of the 
flare and rapidly cooling afterwards. The size of the affected re- 
gion in the photosphere agrees very well with the area derived 
in Sect. 16.11 Our chromospheric modeling in Paper 11 confirmed 
the strong heating of the lower atmospheric layers, during the 
flare the chromosphere is shifted towards higher column masses, 
i.e. deep into the photosphere. In Paper 1 we presented a de- 
tailed list of emission lines observed in the optical flare spec- 
tra; the included lines provide good indicators for the structure 
of chromosphere and transition region. However, as shown in 
Paper II, a self-consistent model of the flare and especially of 
its early phases would require to consider also non-equilibrium 
conditions. Additionally, the flare plasma is non-static, line shifts 
and asymmetries are tracers of mass motions, i.e. chromospheric 
evaporation into the corona and later on material that has cooled 
down and falls back towards the surface. Unfortunately, the lim- 
ited spectral resolution of the X-ray data does not allow to verify 
this for the hot plasma in the corona. 

Nevertheless, our multiwavelength study has revealed a very 
detailed picture of the flare. We were able to determine plasma 
temperatures, densities and abundances, apart from the latter 
not only in the corona but also in photosphere and chromo- 
sphere. We learned that the active region and its associated coro- 
nal loop structure where the flare originated were rather com- 



pact but probably also much more complex than a simple sin- 
gle loop with two footpoints. Multiwavelength observations can 
provi de a comprehe nsive picture of stellar activity phenomena. 
Wo lter et al.l (l2008h have complemented their X-ray observa- 
tions of the fast rotator Speedy Mic with Doppler imaging infor- 
mation, which allowed to localize not only spots on the stellar 
surface but also prominences and a flare. As theoretical models 
become more and more complex, with both stellar atmosphere 
codes and hydrodynamic simulations approaching the non-static 
3D domain, including magnetic fields or further physical condi- 
tions, only detailed observing data can provide complementary 
information. For stellar flares, this requires good spectral and 
temporal resolution over a wavelength range as broad as possi- 
ble. 



7. Summary and conclusion 

We have analyzed simultaneous X-ray and optical data cover- 
ing a giant flare on the active M dwarf CN Leo with regard 
to the plasma properties in the corona and the geometry of the 
flaring structure. As it could be expected from other events of 
similar strength observed in X-rays, temperature and density 
are enhanced by more than an order of magnitude compared to 
the quiescent corona at the flare peak. The flare plasma, which 
can be considered to consist mostly of material evaporated from 
chromosphere or photosphere, shows a different composition. 
Despite its high amplitude, the flare was only of short duration, 
indicating a rather compact flaring structure, which is also com- 
firmed by simple loop modeling. Sustained heating during the 
decay as well as substructure in the lightcurve and hardness ra- 
tio however suggest an arcade-like structure much more com- 
plex than a single flaring loop. Broad multiwavelength coverage 
allowed us to characterize this exceptional event in great detail. 
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